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Abstract: Scalar field theory with an asymmetric potential is studied at zero tempera- 
ture and high-temperature for (ft 6 potential. The equations of motion are solved numerically 
to obtain 0(4) spherical symmetric and 0(3) cylindrical symmetric bounce solutions. These 
solutions control the rates for tunneling from the false vacuum to the true vacuum by bubble 
formation. The range of validity of the thin-wall approximation (TWA) is investigated. An 
analytical solution for the bounce is presented, which reproduces the action in the thin-wall 
as well as the thick-wall limits. 
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1 Introduction 

The problem of decay of a metastable state via quantum tunneling has important applications 
in many branches of physics, from condensed matter to particle physics and cosmology. The 
tunneling is not a perturbative effect. In the semi-classical approximation, the decay rate 
per unit volume is given by an expression of the form 



where Se is the Euclidean action for the bounce: the classical solution of the equation of 
motion with appropriate boundary conditions. The bounce has turning points at the config- 
urations at which the system enters and exits the potential barrier, and analytic continuation 
to Lorentzian time at the exit point gives us the configuration of the system at that point 
and its subsequent evolution. The solution of the equation of motion looks like a bubble in 
four dimensional Euclidean space with radius R and thickness proportional to the coefficient 
of the symmetry breaking term in the potential. When there are more than one solution 
satisfying the boundary conditions, the one with the lowest Se dominates equation (1). The 
prefactor A comes from Gaussian functional integration over small fluctuations around the 
bounce. The zero-temperature formalism is well-developed [1, 2, 3]. In particular, it has 
been proved rigorously that the least action is given by the bounce which is 0(4) invariant 



Linde [4] extended the formalism to finite temperatures. He suggested that at tempera- 
tures much smaller than the inverse radius of the bubble at zero-temperature, the bounces 
are periodic in the Euclidean time r direction and widely separated. Beyond this temper- 
ature they start merging into one another producing what is known as "wiggly cylinder" 

*E-mail : hwidyan@ictp.it, widyan@ahu.edu.jo 



r 





[3]. 



1 



solutions. As one keeps increasing the temperature these wiggles smoothly straighten out, 
and the solution goes into an 0(3) invariant cylinder (independent of Euclidean time r) 
solution that dominates the thermal activation regime. 

A numerical and analytical calculations of the first and second order phase transitions has 
been considered by many authors. For example, an analytical calculation of the nucleation 
rate for first order phase transitions beyond the TWA for the standard Ginzburg-Landau 
potential with asymmetric term has been studied by Miinster and Rotsch [5]. We have 
considered in an earlier work the 4 theory with different symmetry breaking terms [6], 
where we have obtained numerical as well as analytical solution for different values of the 
asymmetric term. In this paper we consider 6 potential motivated by the recent work on 
baryon asymmetry in the standard model with a low cut-off [7]. Also, if the Higgs potential 
is stabilized by a 6 interaction, a strong first order transition can occur for Higgs masses 
well above 100 Gev [8, 9, 10]. Moreover, the 6 potential has been investigated by many 
authors in the context of condensed matter as well as particle physics (see for example 
[11, 12, 13, 14, 15, 16, 17, 18]). 

The general form of the potential is 



which has a second-order transition in m 2 if A > by ignoring corrections due to fluctuations, 
a first-order transition in m 2 if A < 0, and a tricritical point at A = [15]. Since we are 
interested in the case of getting bounce solution, so we take the case of A < 0. Following 
[11], we rewrite the potential U(<f>) in terms of the parameters O and 8 such that 



where 0q = —X/2g and 8 = {X 2 /g — 2m 2 )/4. Looking carefully to the potential, we realize 
that by fixing A and g, then 8 is changed by changing the value of m 2 . Hence the 8 term 
plays the asymmetric part of the potential and is responsible for the first-order feature of 
the phase transition, by causing the coexistence of two minima (false and true) separated by 
a barrier. So, for different values of 8, we get different shape of the potential. 

An interesting special case is the so called thin- wall approximation (TWA), when the 
bubble radius R is much larger than the thickness of the bubble wall and the barrier between 
the two minima is large. In this limit (8 — > 0), there is an analytical formula for Se in terms 
of the wall surface energy, and the details of the field theory are unimportant. However, it 
would be nice to also have an analytical interpolating form for the solution itself. Also, it 
is not clear a priori what the limit of validity of the TWA is. Another interesting case is 
called the thick wall which is reached when the barrier is small. One can easily show that the 
barrier is completely disappeared when 8 = g0 4 and in this case there is no bubbles formed 
and the field goes from the false vacuum to the true vacuum without tunneling. 

In this paper we address the above issues. We obtain accurate numerical solutions for the 
zero-temperature and high-temperature bounces for 6 theory with <p 2 symmetry-breaking 
term. We compute the actions in each case, and find that, for a modest value of the asym- 
metric coupling 8(= 0.1), the action given by TWA formula agrees to within 12.8% with 
that obtained from the numerical solution. We test the criterion for the goodness of TWA, 
in terms of the temperature Tp at which the actions of the 0(4) and 0(3) solutions become 
equal [6]. A numerical investigation shows that the TWA holds up to 8 ~ 0.25. Finally, 
we present an analytical solution which satisfies the equation of motion with parameters 



U(</>) = l -m 2 <p 2 + A0 4 + g<f>*, 



(2) 



W)=^0 2 (0 2 -^) 2 -50 2 , 



(3) 
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fixed by demanding stationary action. This reproduces TWA results very well and, in the 
thick- wall limit, is in good agreement with the numerical results. 



2 Bubble formation 

Let us consider a scalar field theory with a Lagrangian density 



C{<t>) = -{d,<t>) 2 -U{<t>), 



(4) 



where the potential U((/)) has two minima at 0_ (false vacuum) and + (true vacuum). 

In the semi-classical approximation the barrier tunneling leads to the appearance of 
bubbles of a new phase with = <p + as classical solutions in Euclidean space (i.e., imaginary 
time r). To calculate the probability of such a process in quantum field theory at zero 
temperature, one should first solve the Euclidean equation of motion : 



dU{<p) 
dcf) ' 



(5) 



with the boundary condition <fi — > 0_ as x 2 + r 2 — > oo , where r is the imaginary time. The 
probability of tunneling per unit time per unit volume is given by 

o-Se[4>] 



T = A 



(6) 



where Se[4>] is the Euclidean action corresponding to the solution of equation (5) and given 
by the following expression : 



d 4 x 



(7) 



It is sufficient to restrict ourselves to the 0(4) symmetric solution <j)(x 
is this solution that provides the minimum of the action Se 
takes the simpler form 

d 2 (f) t 3d(f) dU{<j>) 



since it 



dp 2 ^ p dp 

where p = ^/x 2 + r 2 , with boundary conditions 



d<j> 



[3]. In this case equation (5) 

(8) 



as p — > oo, 



-f = at 
dp 



p = 0. 



(9) 



We denote the action of this solution by S4. 

Now let us consider the finite temperature case. Following [4], in the calculation of the 
action Se(4>) the integration over r is reduced simply to multiplication by T -1 , i.e., Se[4>] = 
T _1 £3 [0]. Here Se[4>] is the four- dimensional action and S 3 [(p] is the three-dimensional action 
corresponding to the 0(3)-symmetric bubble and given by : 



,[*] = / 



d 3 r 



-(V0) 2 + f/(0,T) 



To calculate Ss((f)) it is necessary to solve the equation 

d 2 (j> 2 d(f) _ dU(<j>, T) 
dr 2 r dr dcf) 



(10) 



(11) 
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with boundary conditions 



d<f> 

<p — > 0_ as r — > oo, — = at r 

ar 



0. 



(12) 



where r = The complete expression for the probability of tunneling per unit time per 
unit volume in the high-temperature limit (T » R^ 1 ) is obtained in analogy to the one 
used in [2] and is given by: 

T(T) = A(T) e~ s ^/ T . (13) 

In the theory of bubble formation , the interesting quantity to calculate is the probability of 
decay between = <p_ and <fi = <p + which are the two minima of U(<f>). There is an interesting 
case (in the sense that the action can be calculated analytically) when U (<fi+) — U ((f)-) = e is 
much smaller than the height of the barrier. This is known as the thin-wall approximation 
(TWA). At T = 0, in the TWA limit, the action S4 of the 0(4)-symmetric bubble is equal 
to 



2£ 3 ■ 

Here Si is the bubble wall surface energy (surface tension), given by 



POO 

Si = / dp 
Jo 



dp' 



(14) 



(15) 



and the integral should be calculated in the limit e — > 0. The bubble radius R is written in 
terms of Si and e as 

R=^. (16) 

£ 

The results presented above were obtained by Coleman [2]. 

These results can be easily extended to the case T » R^ 1 [4]. To this end it is sufficient 
to take into account that 



roo 

An / dr r 2 
Jo 



2 K dr 



) 2 + W,T) 



r e7r J R(T) 3 + 47T J R(T) 2 Si(T), 



(17) 



where Si(T) is the bubble wall surface energy (surface tension) at finite temperature and is 
given by: 



Si(T) 



dr 



<P 2 + W,r) 



As before, the integral should be calculated in the limit e — > 0. 

The bubble radius R(T) is calculated by minimizing S3 with respect to R(T) and this 
gives us 

2S X (T) 



R(T) 



whence it follows that 



lQnSfjT) 
3? ' 



(19) 
(20) 
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3 Numerical results 



For 0(4) symmetry at T = 0, equation (7) reduces to 

"1 



/■CO 

= 2tt 2 / dp , 
Jo 



2 K dp 



r + u{<t>) 



(21) 



We compute the action for different values of the parameter 8 in the symmetry-breaking 
term in the potential £7(0), equation (3), which reads as 



S 4 = 2ir 2 / dp p 6 
Jo 



- 2 (f p ) 2 + 9<P 2 ^-^o) 2 -H 2 



(22) 



Following [11], we assume O — 2.39 and g = 0.07, then the only adjustable parameter in 
the Lagrangian is 8. So, by covering the whole range < 8 < #0q we should be covering all 
relevant cases. 

The equation of motion is now 

d\ 



+ 



3 d<fi 
dp 2 ' p dp 

and the boundary conditions are 



= as p — > oo, 



dp 



at p = 0. 



(23) 



(24) 



By solving equation (23) numerically for different values of 8, substituting the solution in 
equation (22) and integrating, we obtain the action for each value of 8. 

At high temperature, we look for the 0(3) symmetric solution with cylindrical symmetry. 
Then equation (10) takes the form 



S, = 4tt 



f 

Jo 



dr r^ 



J) 2 + 9<P 2 (<P 2 -4) 2 -S<P 2 



The equation of motion is then 



dr 2 r dr 



8^ 3 + 2(# 4_ 5)0) 



and the boundary conditions are 

= 



as 



oo, 



dr 



at r = 0. 



(25) 



(26) 



(27) 



Again, we solve equation (26) numerically for different values of 5, substitute the solution 
in equation (25) and integrate to obtain the action for each 5. Figure 1 shows the bubble 
profile for different values of 8. Note that the value of the scalar field inside the bubble 
decreases with 8. In figure 2 we have plotted this value together with the minimum of the 
potential U((f>). At 5 = the value of 0o, i.e. 0(r = 0), coincides with the minimum of the 
potential m . However, as 8 decreases, the minimum increases while O initially increased 
then it decreases and moves away from the minimum of £7(0). Same behavior has been 
obtained and explained by [19] and it is due the decreasing of the height of the potential and 
the increasing in the energy difference between minima. So, physically this means that as 
the barrier between minima disappears, it becomes easier to from a large bubble with a small 
value of inside it. Same result can been obtained also for the case of zero temperature. 
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Table 1. Numerical values of the action T = and high temperature for different values of 

the asymmetry parameter 5. 



5 


64 (Numerical) 


S 3 (Numerical) 


0.1 


70978.1 


1620.08 


0.2 


9739.27 


441.89 


0.29 


3625.65 


225.41 


0.4 


1519.96 


127.73 


0.6 


523.12 


63.36 


0.8 


253.85 


38.72 


1.0 


143.42 


26.27 


1.2 


r\ -1 r\>-7 

91.07 


18.98 


1.4 


60.85 


14.08 


1.6 


42.67 


10.43 


1.8 


30.03 


7.79 


2.0 


21.34 


5.53 


2.2 


15.01 


3.52 


2.28 


12.84 


2.89 



As we discussed in the introduction, for small values of 5 we can use the TWA formula 
for computing the action. From equation (15) 



S\ = I dr 



DC 







{-^) +g<t>{<t> -<t> Q ) 



= 3.05, (2* 



for O = 2.39 and g = 0.07, see [11]. The radius is given by 

R= 3 -^, (29) 

where e = 4>18 (see [2, 11]). For 5 = 0.1, we have R = 16.02 and the value of the action is 
(see equation (14)) 

S 4 = 61918.4. (30) 

Comparing this analytical value with the numerical one for 5 = 0.1, we get an error equal 
12.8%. In [11], the authors choose 5 = 0.29 to represents the TWA and they have concluded 
that it is not a good value to be taken. This is also confirmed by our calculations where we 
get an error approximately 30%. 

At high temperature, again Si = 3.05. The value of R(T) = 10.68 at 5 = 0.1 and the 
action is (see equation (20)) 

S 3 (T) = 1454.26. (31) 

Comparing this analytical value with the numerical one for 5 = 0.1, we get an error equal 
10.24%. Thus even for S as small as 0.1 the TWA formula for the action does not give very 
accurate results. Obviously, there is no point in comparing numerical results obtained for 
higher values of 5 with the TWA formula. 
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To test our numerical method (we have used Hamming's modified predictor-corrector 
method for solving the equation of motion), we have calculated the action S4 for small 
values of the symmetry breaking parameter 8 in the potential and compared it with the 
TWA formula. In figure 3, we plot the percentage error in the TWA formula as a function 
of 5. The crosses represent our results while the solid line shows a fit to the data. We see 
that the error decreases for small 5, as expected, and approaches zero as 5 — > 0. 

As already mentioned, at zero temperature the 0(4) symmetric solution has the lowest 
value of Se, i-e., Se = S4. At high temperature, we have Se = S3/T. At intermediate tem- 
peratures other solutions exist. In the TWA, however, it has been shown [20] that all other 
solutions have higher Euclidean action. This corresponds to a first order phase transition 
from quantum tunneling at low temperature to thermal hopping at high temperatures. The 
transition temperature Tp is given by equating S4 with S 3 /T, i.e., 

Tp = ^ (32) 

If the surface tension Si is temperature independent, we have 

_ 27^ 



(33) 



Dividing equation (33) by equation (34) and putting e = 4>1 5 (see [2]) we get 

Tp = C*5 (35) 

where 

Thus we see that, in the TWA, Tp increases linearly with 5. We test this by computing 
S3/S4 from our numerical solutions at different values of 5. Figure 4 shows our results for 
the potential given by equation (3). We see that, for S < 0.2, there is very good agreement 
with the predicted linear dependence. This also confirms that, in the domain of validity of 
the TWA, the surface tension Si(T) is independent of T. Beyond 8 ~ 0.2 in our dimensionless 
units, there is a systematic deviation from linearity. Thus we can say that, for values of 5 
larger than this, the wall thickness becomes important. Same behavior has been obtained 
also in our earlier work [6]. 



4 Analytic solution for zero temperature 

We calculate the action analytically in two extreme limits: the thin-wall and thick-wall using 
the potential given by equation (3). 
Thin-wall limit : 5 — » 

In an earlier paper [6], we have found that an analytic solution for the bounce of the 
form of a Fermi function: 

^P) = e (p*-i?VA2 + 1 ( 37 ) 
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is a good approximation for the 4 theory. But it has been shown that for the 6 potential, 
the analytic solution for the bounce has the form [13, 14] 

^aii(p) = T^j, (38) 

where \i = \^8g<pl = 4.21, and p? is the second derivative of the potential in the TWA limit 
evaluated at 0o- So, motivated by the above results, we assume 

^(P) = e (pS-#)/A2 +1 ' ( 39 ) 



where p = y/x 2 + r 2 , R is the radius of the bubble and A its width, acts like a bounce in the 
TWA and leads to the correct value for the action S4. The parameter ^7 is approximately 
equal to true minimum in the TWA. The bounce has values = v /7 at p = and at 
p — > 00. The boundary conditions (9) are satisfied by equation (39). 

To evaluate 7, R, and A, we substitute the ansatz (39) in equation (23) : 

+ ~ p f p = Qgf - 8^0 3 + 2(g4 - 5)0. (40) 

Then the left-hand side (L.H.S.) and the right-hand side (R.H.S.) are respectively 
T n o (3p7A 4 )y/7 , (-4p 2 /A 4 + 4/A 2 ) ^7 

jL - jH ' D - ( 6 (P 2 -R 2 )/A 2 + 1)5/2 + ( e (p 2 -R 2 )/A2 + 1 )3/2 

(p 2 /A^-4/A 2 )V7 , 41 , 

^ ( e (p 2 -R 2 )/A 2 + 1)1/2- ^ 



M.S. = %75/2 %0 ° 73/2 



( e (p 2 -i? 2 )/A 2 + 1)5/2 ( e (p2-R2)/A2 + 1 )3/2 

^ ( e (p 2 -i? 2 )/A2 + 1 )l/2- ^ 

In the TWA, the solution is constant except in a narrow region near the wall at p = R. So, 
we replace in equation (41) 

3p2 3R 2 l 

— r by — — -(1 — aA 2 /R 2 ) in the . -. 9 ttttt term, (43) 

A 4 A 4 ( e (p 2 -R 2 )/A 2 _|_ 1)5/2 ' V ) 

4 4p 2 , 4R 2 , a2/ „ 2 , . , 1 
^ - T^T by - — (1 - bA /R ) m the (e(p2 , R2)/A2 + 1)3/2 term, 

p 2 4 R 2 1 
A* ~ A 2 " by A^ (1 - cA V r2 ) in the (e(p2 -R 2 )/A 2 + ^1/2 term > 

where a, 6 and c are parameters to be determined later. 



(44) 
(45) 
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Comparing equation (41) with equation (42) in the range R 2 (l — A 2 /R 2 ) = R 2 — A 2 < 
p 2 < R 2 + A 2 = R 2 (l + A 2 /R 2 ) where p 2 ~ R 2 as A 2 / R 2 « 1 , we have : 

2 91 2 = ^(l-aA 2 /R 2 ), 



2^7 = ^(1-^7^), 



(46) 



2{g<t>*-8) = ^{l-cA 2 /R% 
We can now evaluate the zero-temperature action S4 : 



Sa = 2tt 2 



dp p 3 



2 y dp 



f + U&) 



(47) 



Substituting equation (39) in equation (47) and integrating we get 

r 2 A 2 

V in, ... - d) i 



S 4 = 2vr 2 J R 4 7 
1 



I 1 7T 2 A 2 . ^ fl 7T 2 A 4 \ 



3A 2 4i? 2 24i? 2 



12i? 4 / 



2A 2 tt 2 A 4 * 
+ 



#o7 1 



# 2 3i? 4 J 
37^ tt 2 A 4 * 

^ 2 + ( 3 + ^4 



(48) 



We now determine the parameters a, 6, and c by demanding dS^/dR = dS^/dA = dS^/dj = 
0. Differentiating equation (48) and using equation (46), we find that to leading order in 
A 2 /R 2 , 



3 - 2a + 46- 2c = 0, 
2 + 3a - 46 = 0, 
-3a + 46 - c = 0, 



(49) 



which leads to a — — 1, b — — 1/4 and c = 2. Using equation (46), we can rewrite equation 
(48) as : 



S 4 = 



ir 2 R 6 



g<& A6 
+ 



1 a b c\ I 7 11a 76 

4-8+1-8 +<-+ 



7r 2 a 



3 tt 2 a 
32 + 24 ~ 32 ~ ~24 



6 

~ 8 + 12 



c_\ A^ 

16 ' 32 16 32 J R? 

7T 2 6 7T 2 C\ A 4 



+ I 96 32 96 + 48 nc ' " fi + ( " 8 ^ 



24 y i? 4 

A 8 
- Of— 
96 / R 6 ' l i? 8 



(50) 



This gives 



7r 2 R 6 / A 2 A 4 A 6 A 8 

; •' 0.063 + 0.141-- 0.049- - 0.020- + O(-) 



b 2 A6 



(51) 
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The quantities 7, R and A are determined from equation (46) using the values of a, b, and 
c. So we have 

c — a\ n , , 4 „ s (b — a 



which gives 



2 91 2 - 2^7 ( + 2 W " ^) = °> ( 52 ) 



R 2 _ gcficry - bjgg - 5) 
A 2 <tf§ 7 - (g<f>t - 5) ' 



(53) 



and 



A 2 — ^ — a _ c — 6 

"2^-2^7" 2^ 7 -2(^-5)' 1 j 

with 7 given by equation (52). We have then, for 5 = 0.1, 7 = 5.83, which implies that 
i^/A 2 = 35.1, R = 16.3, A = 2.75 and S 4 = 70997.3. Comparing these results with the 
TWA formulae, we find that the departure of the radius from the TWA is R/Rtwa = 1-02 
while the departure of the action is S^/Stwa = 1-14, which is a fairly good result. On the 
other hand, there is no departure of the radius as well as the action from the numerical 
values at 5 = 0.1 which is an excellent result. Table 2 shows our numerical as well as the 
analytical values of the action and the radius for different values of 5. We have calculated 
the numerical value of the radius when the derivative of the filed is maximum while in [19] 
the author has calculated the radius in a different way. 

Table 2. Numerical and analytical values of the action and the radius for different values 

of 5. 



5 


(Numerical) 


^4 (Analytical) 


R (Numerical) 


R (Analytical) 


0.1 


70978.1 


70997.3 


16.3 


16.3 


0.2 


9739.27 


10008.3 


8.3 


8.28 


0.29 


3625.65 


3622.26 


5.8 


5.79 


0.4 


1519.96 


1540.44 


4.2 


4.26 


0.6 


523.12 


543.96 


2.86 


2.9 


0.8 


253.85 


266.95 


2.2 


2.22 


1.0 


143.42 


156.01 


1.82 


1.81 


1.2 


91.07 


101.54 


1.53 


1.53 


1.4 


60.85 


70.97 


1.32 


1.33 



Notice that there is an excellent agreement between the radii while actions are fairly 
agree till 5 = 1.0. So, we conclude that our ansatz gives us far better results than the TWA 
formula. In figure 5 we compare our numerical result with the analytic one for 5 = 0.1. From 
the figure we see that the Fermi function agrees very well with our numerical results 

Thick-wall limit: 5 — > g<pQ 
The form of the bounce in equation (39) suggests that the thick wall limit, which would 
correspond to small values of R 2 /A 2 , would be obtained by approximating the Fermi function 
by the Maxwell-Boltzmann function, which leads to a Gaussian: 

2 (p) = ie - p2/A2 . (55) 
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The action for this form of bounce is found to be 

5 4 = 2tt 2 A 4 7 



(56) 



Equations (46) then reduce to 

a 



(57) 



Note that in this case 7 <C 1, so 7 2 is negligible (a = 0). 

The values of b and c are again obtained by demanding dS^/dA = dS±/d^ = 0. The 
relation between them is 

2 + 6- c = 0, 
2 + b - 2c = 0. 



This gives b = —2, c = 0, giving 



This yields the action 



5 = #0, #o7 



25 



A 2 ' 

'•' + 0(--J = 12.37 



(58) 



(59) 



for 5 = 2.28. The numerical value is S 4 = 12.84, so the error is 4.6%. 

Thus, the form of the bounce given by equation (46) seems valid over the whole range of 
5 (from to 2.28), and in the two extreme limits is amenable to analytic calculations. 



5 Analytic solution for high temperature 

We discuss now the high-temperature action S3 for the thin wall limit as well as thick wall. 
Thin-wall limit : 5 — » 

The bounce takes the following from: 



e (r 2 -R2)/ A 2 _|_ j ) 



(60) 



where r 2 = x 2 and the other parameters R and A have the same physical significance in three 
dimensions. The boundary conditions given by equation (27) are satisfied by the bounce. 

We substitute the bounce in the equation of motion (26) and assume the solution is 
constant except in a narrow region near the wall r = R. The resulting equations enable us 
to evaluate the action given by equation (25), and after integrating we get the following: 



47Ti? 3 7 



8A 2 16R 2 



a2 (2 A 2 tt 2 A 4 \ 2 /l 
- #o7 3"^ + l2^ + 91 3 



3A 2 tt 2 1 a 4 ; 

4R 2 + ^24 + 8^ 4 



(61) 
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In terms of parameters a, 6 and c, the action takes the simpler form 

4tt R 5 



2g<j% A 6 



1 a b c 
4 ~~ 6 + 3 ~~ 6 



/9 5a 26 7c\ A 2 
+ \8~24 + y~12j^2 



/7 7T 2 5a 
+ 8 + 64 + T 



7r 2 a 



48 



7b 7T 2 b 

T + 1m 



+ 



7tt 2 



7a 
32 



7vr 2 a 7tt 2 6 
+ 




7n 2 c\ A 6 ^A 8 , 



^128 32 96 128 96 / 
where the relations between a, 6, and c to leading order in A 2 /R 2 are 

1 -a + 26-c = 0, 
1 + 3a - 46 = 0, 
-3a + 46 + c = 0, 



(62) 



(63) 



which leads to a = —5, 6 = —7/2 and c = —1. Hence the action in equation (62) is reduced 
to 

4-7T R 5 ( A 2 A 4 A 6 A 8 \ 

s ' = m f l - 083 + °- 417 ^ + °' 699 ^ + °- 914 ^ + 0( ^ ) ) • (64) 

where 7, A and R are obtained from equations (52), (53) and (54). 

We have then, for 5 = 0.1, 7 = 6.04, which implies that i? 2 /A 2 = 22.62, R = 11.06, 
A = 2.33 and S3 = 1691.95. Comparing these results with the TWA formulae, we find that 
the departure of the radius from the TWA is R/Rtwa — 1-02 while the departure of the 
action is S^/Stw — 116, which is a fairly good result. On the other hand, there is a very 
small departure of the radius as well as the action from the numerical values at 5 = 0.1 
which is an excellent result. Similarly as in the case of zero temperature, if you go to higher 
values of S, then there will be a departure from the numerical results. 

In the TWA the radius of the bubble is much greater than its thickness. So, for 5 = 0.1, 
we get A = 2.75 which is much less than R = 16.3 as is expected. Same result is obtained 
for the zero temperature as well as in [19]. 

Thick-wall limit : 5 — > g<pQ 

At higher temperature the bounce takes the form as the case of zero temperature, 

2 (r) = 7e - r2/A2 , (65) 

with the action 



S 3 = 7rA 3 v / 7r7 



3 1 4 
Ta^ + - ^0o7 + (#0 - 5) 



(66) 



Again defining 1g<$fi = — 6/A 2 , 2(g0Q — 5) = — c/A 2 and neglecting 7 2 , we find 6 and c as 
before by demanding dS^/dA = dS^/dj = 0. The relation between 6 and c is given by 



3 

16 



2 5 /2 8 



3c 



3 36 

Tq + ¥J 2 ~~8 



0, 



0, 
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which leads to b = —y/2 and c = -1/2, giving A 2 = 1/(4^ - 45) and 7 = (4^0^ - 
A5)/{y/2g<j>l). The action can be simplified to 



Note that the value of the action is independent of A and depends only on S and if 5 = g<pQ, 
then 7 = and the action S3 = which is consistent with the result that the hump of the 
potential will disappear at this value of S. Another important result is that A will diverge 
in the limit S = g$\ which has been also obtained by [19]. So, to get a real value of action 



We have plotted in figure 6 the numerical and analytical bubbles for S = 2.0. Note that 
in spite of the discrepancy in the value of 7 for the numerical and analytical profiles which is 
due to the neglecting the terms of order R 2 /A 2 in equation (67) the departure of the actions 
is small, i.e. S3 (numerical) / S3 (analytical) = 1.05. 

6 Conclusions 

We have obtained accurate numerical solutions for the zero-temperature and high-temperature 
bounces for 6 potential with 2 symmetry-breaking for the entire wall thickness interval 
< 5 < g<pQ- We compute the actions in each case and find that, for a modest value of the 
asymmetric coupling 5 = 0.2, the action given by the TWA formula agrees to within 12.8% 
with that obtained from the numerical solution. At high temperatures, the conclusion is 
qualitatively similar. 

We have checked our numerical method by comparing the action obtained numerically 
with the one obtained from the TWA formula. Very good agreement is obtained as we go to 
small values of 5. We also verify that as 5 is reduced the error in the TWA formula goes to 
zero. We check the criterion for the goodness of TWA proposed in [6], in terms of the relation 
between S and the temperature Tp at which the actions of the 0(4) and 0(3) solutions become 
equal. A numerical investigation shows that TWA holds up to 5 ~ 0.2. Finally, we present 
an analytical solution which satisfies the equation of motion in an approximate sense in two 
limiting cases. The first of these reproduces the leading corrections to the TWA results very 
well and it fairly matches the numerical results of the action up to 5 = 1.0. The second is 
applicable for the opposite case of a very thick wall. This gives us insights into the nature 
of the bounce solutions for various values of 5 going from thin to thick walls. 

Some of our results match very well with those obtained in [19]. For example, we get 
the same behavior of the the minimum of the potential <p m and the value of the inside 
the nucleated bubble, O , see figure 2. Moreover, the divergence of the thick of wall at the 
vanishing of the hump of the wall is obtained in [19] numerically while we get the same 
behavior analytically. 

Much of the work on inflationary models relies on the zero-temperature potential, so 
our results could be relevant for inflation [21]. They may also have some bearing on the 
formation of topological defects in a first order phase transition where authors consider 
zero-temperature potentials, see for example [22]. 

So far, we have discussed the action only at zero and high temperatures. To obtain the 
bounce solution at intermediate temperatures, we have to solve a partial differential equation 




(67) 



we must always have 5 < g<pQ. 
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with periodic boundary conditions in the r direction either numerically or analytically. This 
work will be presented in a future publication. 
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Figure Caption 



Figure 1. Shape of the critical bubble at different 5. 

Figure 2. The minimum of the potential 4> m and the value of the inside the nucleated 
bubble, 0o- 

Figure 3. Error in the TWA formula as a function of 5. The crosses represent our results 
while the solid line shows a fit to the data. 

Figure 4. Deviation of Tp from the TWA limit. The dashed line represents the TWA 
limit while the crosses are our numerical results. 

Figure 5. as a function of p. The dashed line is the Fermi function while the doted 
line is the numerical result. 

Figure 6. as a function of r. The dashed line is the Gaussian function while the solid 
line is the numerical result. 
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